Lefschetz-thimble & complex Langevin approach 
to Silver Blaze of one-site Hubbard model 


Yuya Tanizaki 


Department of Physics, The University of Tokyo 


Theoretical Research Division, Nishina Center, RIKEN 


Dec 3, 2015 O KEK 


Collaborators: Yoshimasa Hidaka, Tomoya Hayata (RIKEN) 


a e = = Ss nac 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Dec 3, 2015 @ KEK 1/18 


Finite-density QCD? 


Neutron star 
e Cold and dense nuclear matters 
o 2Msun neutron star (2010) 
o Gravitational-wave observations (2017) 


Reliable theoretical approach to equation of state must be developed! 


Z(T,u) = f DA Det(P(A, u) +m) exp —Sym(A). 


quark gluon 
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Silver Blaze problem in finite-density QCD 
QCD & |QCD| 


Zoco = [PA (dety,D,) gr P. Z\acp! = [ra |dety,D,| B. 


At u = 0, these two are the same! But, 


ng 
QCD 


OCH 


m /2 my/3 "H 


|QCD] experiences the phase trans. at y = m,/2 ~ 70MeV. 
In QCD, the state = QCD vacuum for u < my /3 ~ 300MeV. 
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Complexification of fields 


There are two “new” approaches to the sign problem: 
e Complex Langevin method: Solve the Langevin eq. 


1] is a real stochastic noise, (n(0)n(0')) = 20(0 — 6’). 
We have no solid foundations. 


e Path integral on Lefschetz thimbles: 


f din el) — Sitka) f digg P 


oO 


Path integral is performed on steepest descent paths JJ 
This is mathematically rigorous. 
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Lefschetz-thimble integrations 


Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = | evi E + 2 


Complexify the integration variable: z = x + iy. 


Im(6) 


4 


Integrand on R, and on A 
(a=1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
d 3 
Ai(a) = > No [ Pi E taz) 


ng: intersection number of the steepest ascent contour Ke and R. 


Figure: Lefschetz thimbles 7 and duals K (a = 100%, =i) 
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Sign problem & One-site Hubbard model 


One-site Fermi Hubbard model 
One-site Hubbard model: 
A = Um, — u(fu + fy). 
Fock state gives the number density immediately: 


10 Zielt + ePQu-U)) 
bôu id= 1 + Zen + ed2u-U) 


In the zero-temperature limit, 
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Sign problem & One-site Hubbard model 


Path integral for one-site model 


The path-integral expression for the one-site Hubbard model: : 


B d A 2. UD 
a=a > | ay (14 ef Ge)? 9807/20, 
27U R EE 


Fermion Det 


Integrand has complex phases causing the sign problem. 


( is an auxiliary field for the number density: 
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Sign problem & One-site Hubbard model 


444 
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144 


Flows at y/U < —0.5 and u/U > 1/5 


2 


Re(z) 


Figure: Flow at u/U 


/ dz e°@), 
Se 


Number density: n, = 0 for u/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 


Z- 


9 / 18 


Dec 3, 2015 0 KEK 


Lefschetz-thimble path integral 
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Sign problem & One-site Hubbard model 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at u/U = 0 


py p 1. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


m 


Complex saddle points lie on Im(z,,)/U 
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Lefschetz-thimble path integral 
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CL distribution at —0.5 < u/U < 1.5 


BU =30 u/U=0 


o 


o 
bo 


0.4 


B[Im(z) — (u + U/2)] 


0 
BRe(z) 
Figure: CL distribution at y/U = 0 


It looks quite similar to Lefschetz thimbles. 
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Sign problem & One-site Hubbard model 


Curious incident of n in the one-site model 


We have a big difference bet. the exact result and naive expectation: 


25 OÖ 05 1 15 pfu 


CL reproduces the naive onel(YT, Hidaka, Hayata, arXiv:1509.07146, 1511.02437) 


This is similar to what happens for QCD and |QCD|. 
p/U = —0.5 S u = mz /2. 
(cf. Monte Carlo with one-thimble approx. gives the naively expected results. 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258) 
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Semiclassical analysis 


Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


WU x 1/2 


U 
Ea cR By +2rmT. 


At these solutions, the classical actions become 


Sox -£ (+, 
Re (Sin — Sp) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 
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Semiclassical analysis 


Semiclassical partition function 
Using complex classical solutions Zm, let us calculate 


oo 
La t= d ee, 


M=—-09 


This expression is valid for —1/2 S u/U $ 3/2. 
This is calculable using the elliptic theta function: 


Za 


K 


= 1 
-S [149 9 Heo \ ntm?/pu 
e + 2 am (5 + 2 e 


— eo pe 1 —2n? /BU 
e ^ ((5*5).« j 
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Semiclassical analysis 


Number density & Lee-Yang zeros 
Lee-Yang zeros of Za: 


e e 
e e 
a e Re 
. e 


Semiclassical study gives the correct transition! 


" 2 (1<p/U « 3/2), 
Nel = 39, Za + 1 (0 « p/U <1), 
Pap 0 (—1/2 < u/U <0). 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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Semiclassical analysis 


Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za 7 $e A e 0727). 


o Lee-Yang zeros cannot appear at u/U — 0, 1. 
e One-thimble, or "phase-quenched" , result: n ~ u/U + 1/2. 
e Complex Langevin ~ the phase-quenched multi-thimble approx. 


Consequence 


In order to describe the step functions, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schäfer, Ünsal, PRL 115 041601) 
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Semiclassical analysis 


Complex Langevin simulation 
One-site Fermi Hubbard model: 


=== exact 
2.0} HEH w/o phase 


‚9: w/ phase 

= 1.5 
o 
io — — | E . 

0.5 

0.0 

-0.5 0.0 0.5 10 Ls 
p/U 


(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 


Consequence 


Modified complex Langevin is not perfect yet, but it seems to point a 
correct way. | hope this can attack the Silver Blaze phenomenon. 
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Summary & Conclusion 


Summary and Conclusion 


e Picard-Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


e One-site Hubbard model is a nice toy model to play with the 
sign problem. 


e Destructive and constructive interference of complex phases 
among Lefschetz thimbles play a pivotal role for the baryon 
Silver Blaze. 
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Backups 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles 7,: (classical eom S'(2,) = 0) 


/ d"z e 3) = Y (Ke; y) d'zg 9i 


o 


Js are called Lefschetz thimbles, and Im[S] is constant on it: 


M =}, det _ eo 


1—-—oo dt pe 


a [40 


(Ko, R}: intersection numbers of duals K, and R” 


(Ko = {2(0)|2(00) = 24 )). 


[Witten, arXiv:1001.2933, 1009.6032] 
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Numerical results 


Results for BU = 30: (1, 3, 5-thimble approx.: Jo, Jo U J: 


+1, and Jo U J41 U Jae ) 


2.0 
— exact 
— 1-thimble 


15 --- 3-thimble 


n 
1 
11 
d 
n 
14 
14 
UI 
14 
my 
my 
ry 
na 
A 
1 
1 


(YT, Hidaka, Hayata, arXiv:1509.07146[hep-th]) 
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———————— SE ve: 
Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz,(0) OS 


ag 7 gz 00) + V(O). 


0: Stochastic time, n: Random force satisfying 
Lu Ia Uu = 20(0 0). 
Ito calculus shows that 


LOO), = A(O" (z,(8))), — (O'(Zn(9))S"(2n(9)))n- 


If the I.h.s becomes zero as O > oc, this is nothing but the 
Dyson-Schwinger eq. 
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Complex Langevin and Lefschetz thimbles 


For any solutions of the DS eq, 


(05) 7 7 374. | Mo, 


c o 


ind, € C. To reproduce physics, d, = (Ko, R). 


So far, we ONLY assume the convergence of the complex Langevin 
method. 
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Semiclassical limit 
Let us take h < 1 for computing 


(Ote) = ds [ dzeS/*O(z). 


| have NO idea how to compute the LHS. However, positivity of the 
probability density and its localization around zo's imply that 


"ce. 3 (0%) S y O. 


RHS is 


[ die SOM"O(z) > Sue Oe): 


Yuya Tanizaki (University of Tokyo, RIKEN) Lefschetz-thimble path integral Dec 3, 2015 @ KEK 24 / 18 


Semiclassical inconsistency 
In the semiclassical analysis, one obtains (for dominant saddle points) 


C, > 0. And, d, = (K,,R) to get physics. = Inconsistent! 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points. 
e Those saddle points have different complex phases. 


Consequence 


Naive complex Langevin method cannot explain the Silver Blaze 
phenomenon. 
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Proposal for modification 
Assume that 


_ Ko, R) Zoch So) 


S" (Ze) 


Because of the localization of probability distribution P, it would be 
given as 


P= re supp(P,)Nsupp(P;) = 0. 


CO 


Assumption means “CL = phase quenched multi-thimble approx." : 


(Ko, R) 


27h 
T äich 


Sz) De): 


(Olh =S 


CO 
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Proposal for modification (conti.) 


If so, defining the phase function 


(zc)] -iIm Size 


$(z,z) = Kä GER "Xsupp(P,) (2,2), 
we can compute 


This new one is now consistent within the semiclassical analysis. 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 
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